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ABSTRACT 


Examination  of  the  seismic  source  from  both  the  near-  and  far- 
fields  has  been  undertaken.  Using  a discrete  Fourier  transform  (DFT)  and  a 
maximum  entropy  spectral  estimator  on  near-field  acceleration  data  produced 
corner  frequency  and  low  frequency  level  estimates.  This  latter  estimator 
eliminates  leakage,  thereby  obtaining  more  accurate  values  for  the  high  fre- 
quency enu,  with  better  roll-off  values  and  definition  of  side  lobes.  When  the 
maximum  entropy  spectral  estimator  was  applied  to  data  from  the  Parkfield 
earthquake  and  the  Bear  Valley  event  of  June  22,  1973,  the  lower  frequency 
level  was  lower  than  that  from  the  DFT  in  all  cases. 

Using  spectral  fitting  procedures  on  fundamental  mode  surface 

wave  data  from  two  central  California  events,  generally  close  agreement  was 

obtained  with  source  mechanism  solutions  from  bodywaves.  For  the  Bear 

Valley  earthquake  of  June  22,  1973,  the  seismic  moment  obtained  from  the 

surface  wave  data  was  an  order  of  magnitude  smaller  than  that  obtained  from 

acceleration  data.  In  an  effort  to  reduce  the  scatter  in  M - m plots  for  an 

s b 

earthquake  population,  several  attempts  were  made  to  reduce  tne  variance  in 

M . Both  demultipathing  and  correcting  for  the  radiation  pattern  had  no  effect. 

It  was  concluded  that  the  variance  of  m,  is  the  controlling  factor.  Finally, 

b 

theoretical  higher  mode  surfa.ce  wave  pectra  was  generated  for  a double  cou- 
ple source  in  a layered  half  space.  Both  the  Rayleigh  and  Love  wave  higher 
mode  spectra  were  found  to  vary  shape  dramatically  as  a function  of  source 
depth. 
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SECTION  I 


SPECTRAL  ESTIMATES  OF  NEAR-FIELD  ACCELERATION 


A.  INTRODUCTION 


The  results  of  the  analysis  of  the  Bear  Valley  Earthquake  of 
June  22,  197  3 using  the  Haskell  moving  dislocation  model  (Turnbull  and  Battis, 
1974)  demonstrate  the  high  degree  of  non -uniquenes s of  the  source  parameters 
solution  for  a given  event  using  this  source  representation.  In  this  study  three 
equally  satisfactory  solutions,  judged  solely  on  the  quality  of  the  velocity 
waveform  fit,  were  found.  Two  of  the  solutions  were  essentially  the  same 
except  for  a shift  in  epicentral  location  of  2 kilometers.  The  third  solution, 
however,  was  radically  different  in  both  dip  angle  and  dislocation  amplitude. 
Without  externally  imposed  restrictions  on  the  source  structure,  such  as  fault 
plane  solution,  seismic  moment  and  fault  dimensions,  it  is  obvious  that  fault 
models  generated  with  the  Haskell  model  cannot  be  considered  as  valid  indi- 


cators of  the  true  source  structure  of  an  event. 


While  fault  plane  solutions  cannot  generally  be  obtained  from 
near-field  acceleration  data  alone,  using  techiques  of  spectral  analysis  in 
conjunction  with  a source  model,  such  as  that  given  by  Brune  (Brune,  1970), 
one  can  obtain  estimators  of  the  seismic  moment  and  corner  frequency.  These, 
in  turn,  provide  an  estimate  of  source  dimensions  and  total  dislocation  which 
can  be  used  to  limit  the  range  of  solutions  derived  from  time  domain  wave- 


form modeling  of  the  earthquake  source. 


Therefore,  in  an  effort  to  limit  the  degree  of  non -uniqueness 
of  the  H iskell  source  model,  and  to  analyze  the  consistancy  of  the  source  para- 
meters between  ume  domain  and  frequency  domain  representations  of  the 
earthquake  source,  a program  of  spectral  analysis  of  the  available  near- field 


* 


acceleration  data  was  implimented.  The  initial  results  from  this  work,  using 
the  discrete  Fourier  transform,  were  given  in  a previous  report  (Turnbull 
and  Battis,  1974).  However,  they  have  subsequently  been  found  to  be  in  error. 
This  section  will  discuss  this  error,  and  its  effects  on  the  results,  an'  addi- 
tional techniques  adapted  during  this  report  period  to  provide  better  spectral 
estimation  from  near-field  acceleration  data. 


B.  FAST  FOURIER  TRANSFORM  SPECTRAL  ESTIMATOR 

The  discrete  Fourier  transorm  (DFT)  technique  for  obtaining 
a spectral  estimator  is  well  documented  and  no  discussion  of  the  theory  will 
be  given  here.  It  was  this  form  of  analysis  which  was  initially  programmed 
for  the  spectral  analysis  of  near-field  acceleration  data  (Turnbull  and  B-ttis, 
1974)  In  this  program  the  autocorrelation  function  of  the  data  was  calculated 
as  an  intermediate  step  in  the  process,  with  the  amplitude  spectra  being  eval- 
uated by  Fourier  analysis  of  the  autocorrelation  estimator.  The  rational  for 
utilizing  this  system  was  hat  the  appropriate  filtering  process  could  be  judged 
more  accurately  on  the  basis  of  the  behavior  of  the  autocorrelation  function 
than  directly  from  the  acceleration  data. 

The  use  of  the  autocorrelation  function  to  determine  the  filtering 
characteristics  proved  satisfactory.  However,  a programming  error  associ- 
ated with  the  normalization  of  the  autocorrelation  function  induced  errors  in 
the  evaluated  amplitude  spectra.  These  errors  take  the  form  of  a scale  factor 
dependent  on  the  time  gate  of  the  analyzed  signal,  specifically 

A(f)  = A*(f)  / ( /tTI  ) 

where  A(f)  is  the  correct  amplitude,  A (f)  is  the  original  calculated  ampli- 
tude and  T is  the  signal  duration.  The  effect  is  to  reduce  the  observed 
amplitudes  and  thus  the  low  frequency  level  of  the  spectra.  However,  no 
change  in  the  spectral  shape  occurs. 
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We  have  previously  reported  on  the  source  parameters  of  18 
earthquakes  derived  from  DFT  spectral  analysis  of  acceleromettr  data  (Turn- 
bull,  et  al.  , 1974b).  As  no  change  in  the  spectral  shape  occurs  from  this  er- 
ror, the  evaluated  corner  frequencies  jf^)  are  still  valid.  The  low  frequency 
level  (Q  being  improperly  scaled,  does  change,  in  turn  affecting  the  esti- 
mated seismic  moment  for  each  event. 

The  calculated  seismic  moments  have  been  corrected  and  are 
given  in  Table  1-1.  In  all  cases,  the  seismic  moments  have  been  reduced  and 
generally  conform  better  with  other  estimates. 


C.  MAXIM  JM  ENTROPY  SPECTRAL  ESTIMATOR 

The  use  of  the  DFT  for  spectral  estimation  from  near-field  ac- 
celeration data  has  one  prominent  drawback  effecting  the  spectrnx  ^ valuation 
of  source  characteristics.  Windowing  effects,  in  the  form  of  power  leakage 
from  one  frequency  to  another,  produce  high  levels  of  noise  in  the  DFT  spec- 
trum (Figures  I- 1 , 1-2,  and  1-3).  A common  example  of  this  effect  is  the  spur- 
ious side  lobes  which  occur  in  the  DFT  spectrum  of  a pure  sinusoid.  This  con- 
dition can  be  corrected  by  smoothing  the  spectrum,  but  this,  in  turn,  reduces 
the  prominence  of  any  spectral  peaks.  In  addition,  since  the  smoothing  pro- 
cess can  be  seen  as  a weighted  averaging  process  on  the  spectrum,  the  accur- 
acy of  the  calculated  amplitudes  within  some  localized  bandwidth  can  be  biased 
by  a prominent  peak  or  hole  in  the  spectra  located  near  this  frequency  band  of 
interest. 

Both  effects  of  windowing  are  of  importance  when  one  is  attempt- 
ing to  evaluate  the  spectral  parameters  which  are  considered  characteristic  of 
the  earthquake  source,  specifically  the  corner  frequency,  fQ,  and  low  frequency 

level,  Q . From  previous  work  in  this  area,  it  would  appear  that  the  most 
o 

consistent  method  of  determining  the  corner  frequency  is  to  locate  the  maximum 
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CORRECTED  PARAMETERS  FROM  SPECTRAL  ANALYSIS  OF 
NEAR-FIELD  ACCELEROGRAM  DATA 
(PAGE  1 OF  5) 
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CORRECTED  PARAMETERS  FROM  SPECTRAL  ANALYSIS  OF 
NEAR -FIE  ED  ACCELEROGRAM  DATA 
(PAGE  5 OF  5) 


FREQUENCY  C HZ*  3 
BEAR  VALLEY  5TAT  8 COMP  N45E 


FIGURE  1-1 

NORMALIZED  DFT  ACCELERATION  ENERGY  DENSITY  SPECTRUM 
FOR  N45°E  COMPONENT  OF  STATION  8;  BEAR  VALLEY  EARTHQt  At  F 
OF  JUNE  22,  1973  (HIGH  PASS  FILTER  CORNER  FREQUENCY  AT  0.7  Hz) 
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FREQUENCY  C HZ.  ) 

BEAR  VALLEY  STAT  8 COMP  N45E 


FIGURE  1-2 

NORMALIZED  DFT  VELOCITY  ENERGY  DENSITY  SPECTRUM  FOR  N45°E 
COMPONENT  OF  STATION  8;  BEAR  VALLEY  EARTHQUAKE  OF  JUNE 
22,  1973  (HIGH  PASS  FILTER  CORNER  FREQUENCY  AT  0.7  Hz) 
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NORMALIZED  DFT  DISPLACEMENT  ENERGY  DENSITY  SPECTRUM  FOR 

0N45°E  COMPONENT  OF  STATION  8;  BEAR  VALLEY  EARTHQUAKE  OF 
JUNE  22,  1973  (HIGH  PASS  FILTER  CORNER  FREQUENCY  AT  0.7  Hz) 
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peak  in  the  velocity  spectrum.  However,  if  a smoothed  DFT  spectrum  is  used, 
the  peak  will  be  lowered  in  amplitude  and  its  bandwidth  widened,  making  the 
accurate  determination  of  f more  difficult.  WTiile  not  truely  significant  in 
the  case  of  high  signal  to  noise  data,  this  effect  can  be  important  in  low  sig- 
nal to  noise  cases  such  as  station  2 or  4 for  the  Bear  Valley  earthquake  of 
June  22,  1973. 


Secondly,  it  is  inherent  in  the  methods  of  applying  the  instru- 
ment response  correction  to  obtain  a displacement  spectra  from  acceleration 
data  that  the  spectral  level,  as  zero  frequency  is  approached,  should  go  to 

infinity.  For  an  accelerometer,  the  instrument  correction  to  obtain  displace- 

-2 

ment  spectrum  is  ( 2 7T f ) . If  the  acceleration  spectrum  does  not  approach 

zero  amplitude  at  zero  frequency  in  a manner  producing  a limited  low  fre- 
quency level,  then  at  low  frequency,  the  spectrum  will  be  unbounded.  Win- 
dowing effects  will  then  boost  H and,  if  smoothing  is  attempted  on  the  spectrum, 
a further  amplification  of  the  lo  »'  frequency  level  will  occur.  Depending  on 
the  degree  of  leakage  and  smoothing,  this  may  greatly  effect  the  estimation 

of  seismic  moment  based  on  0 . 

o 

In  an  effort  to  reduce  the  errors  from  leakage  and  smoothing  the 
maximum  entropy  or  Markov  spectral  estimator  was  adapted  for  near-field 
acceleration  data  analysis.  This  method,  developed  by  John  Burg  (Burg,  1967) 
of  Texas  Instruments  Incorporated,  has  the  advantage  over  the  standard  DFT 
that,  in  theory,  it  has  no  windowing  effects  and  thus  no  leakage  (King,  et  al.  , 
1974).  Therefore,  a smooth  amplitude  spectrum  is  obtained  while  maintaining 
high  resolution.  Of  course,  the  accuracy  of  the  spectral  estimation  is  limited 
by  the  degree  to  which  the  processed  signals  match  the  hypothesis  of  the  techni- 
que. 

A simplified  explanation  of  the  maximum  entropy  spectral  esti- 
mator is  given  below.  More  detailed  examinations  of  the  method  and  its  com- 
parison with  other  spectral  analysis  methods  can  be  found  in  papers  by 
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Barnard  (1975).  King,  rial.  (1974).  Ulrych  (1972),  Burg  (1968),  and  Lacoss 
(1971).  The  following  is  a summary  of  Section  III  of  this  report  by  King, 

et  al. 

The  two  basic  equations  of  discrete  power  spectral  evaluation 
(the  square  of  the  amplitude  spectrum)  are  those  which  define  the  discrete 
autocorrelation  function,  «M  r ).  and  the  discrete  power  spectrum,  P(f).  These 

are  given  as 


1 


7 ) 2T  + 


T E 


XNXN+T 


( - oo  < T <») 


N=  -T 


f 


i 


and 

00 

P(f)  = _L  'JT  «f,(r)  cos  (27Tf  rAt),  (Olf  < W = 1/ 2 At) 

r r -oo 

where  At  is  the  sampling  period  of  the  discrete  time  scries.  As  cf>(  r)  is 
known  exactly  for  all  lags,  then  the  power  spectrum  has  infinite  resolution 

within  the  band  limit  W. 

This  requires  that  <f)(  r)  be  known  exactly  for  all  lags.  As  we 
are  dealing  with  a finite  data  set,  this  does  not  hold.  Any  physically  realistic 
values  of  may  be  used  to  supplement  the  known  values,  but  a criteria  for 

judging  what  is  a realistic  value  may  be  hard  to  determine.  In  many  cases 
the  data,  X , beyond  the  known  signal  is  defined  as  zero,  producing  a tapered 
autocorrelation  function  which  introduces  windowing  effects.  Depending  on  the 
transitory  nature  of  the  signal  this  may  or  may  not  be  realistic. 

The  maximum  entropy  estimator  takes  a different  approach 
and  attempts  to  supplement  the  autocorrelation  function  in  some  "optimum" 
manner.  This  optimum  manner  chosen  in  the  maximum  entropy  technique 
is  that  which  maximizes  the  degree  of  disorder,  or  entropy,  of  the  resulting 

power  spectrum.  Given  N known  lags  of  the  autocorrelation  function,  the 
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manner  in  which  the  N + 1 lag  is  evaluated  while  guaranteeing  a non-negative 
definite  autocorrelation  ( motion,  and  without  modifying  the  known  lags,  turns 
out  to  be  the  same  as  the  design  of  an  N + 1 point  prediction  error  filter  |a|. 


•/'(0)  «Mi)  'A(N)  1 


f/'(i)  «/»( o) 


The  spectrum  P(f),  the  maximum  entropy  spectrum,  is  then 


2A t PE 


P(f)  = 


E v 

j=o 


-2i7TfjAt 


Consider  the  optimum  N point  prediction  filter  operating  on 

A A 

the  time  series  X producing  an  output  X . The  difference  in  X and  X is 
IN  N N N 

which  is  the  output  of  the  N + 1 point  prediction  error  filter.  Since  the 

prediction  filter  is  optimum  in  a least  mean  square  sense,  the  expected  value 

of  the  product  of  € with  any  previous  X is  zero. 

N N 


E[  «NX.]  = 0 (i  = 1 , N) 


Suppose  that  the  N point  prediction  filter  was  able  to  predict 
all  the  cor relationed  components  of  X^  If  this  were  true,  then  the  N point 
filter  would  perform  as  well  as  an  infinitely  long  prediction  filter.  The  output 
of  the  N point  filter  applied  to  the  original  data,  X,  would  be  pure  white  un- 
correlated random  noise  with  a power  density  of  2ArPE  . Thus,  all  spectral 
information  of  the  original  time  series  is  contained  in  the  prediction  filter, 
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and  in  effect,  the  autocorrelation  function  is  extended  to  infinite  lags.  1 j 


recover  the  spectrum  of  the  input  time  series,  one  need  only  divide  the  output 
power  spectrum  by  the  power  response  of  the  filter 


It  should  be  noted  that  the  power  spectrum  of  the  input  time 
series  is  for  the  correlated  part  of  the  signal  only.  Thus,  if  one  assumes 
that  the  uncorrelated  part  of  the  signal  is  due  totally  to  noise,  the  spectrum 
represents  the  true  signal. 


In  Figures  1-4  to  1-6  an  example  of  the  maximum  entropy 
spectral  estimates  of  the  Bear  Valley  earthquake  signal  recorded  by  the  N45 
component  at  station  8 is  shown.  For  comparison,  Figures  1-1  to  I- 1 show 
the  spectra  obtained  for  same  signal  using  the  Dl<  T.  In  both  cases,  a high 
pass  filter  with  a corner  frequency  of  0.6  Hz  has  been  applied  to  the  accelera- 
tion data.  In  addition,  the  maximum  entropy  displacement  spectrum,  calculated 
with  the  unfiltered  data  is  given  in  Figure  1-7. 


It  can  readily  be  seen  that  the  maximum  entropy  spectral  esti- 
mator and  the  DFT  are  in  general  agreement  on  the  spectral  shape  and  ampli- 
tudes. As  expected  though,  the  maximum  entropy  spectra  is  considerably 
smoother.  It  appears,  in  a comparison  of  Figure  1-6  and  Figure  1-7  that 
low  frequency  level  has  a drop  in  amplitude  just  before  the  corner  frequency 


as  would  be  exnected  from  Archambeau's  source  model.  However,  more 


analysis  of  other  events  must  be  made  before  this  can  be  confirmed. 


Due  to  the  elimination  of  leakage  by  use  of  the  maximum  en- 


tropy spectral  estimator,  two  addition  features  of  the  spectral  shape  are  more 
reliably  estimated.  First,  in  DFT,  leakage  will  reduce  the  rate  of  high  fre- 
quency roll-off.  Without  leakage,  a better  estimate  of  this  value  is  attained. 
Secondly,  the  high  frequency  end  of  the  spectrum  is  smoother  and  side  lobes 


can  be  detected.  No  work  with  either  of  these  features  has  been  conducted 


for  this  report,  and  it  has  not  been  demonstrated  that  the  side  lobes  are  really 
those  predicted  by  spectral  source  representations. 
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FREQUENCY  C HZ#  I 
STATION  8 COMP  NA5E 


FIGURE  1-4 

NORMALIZED  MAXIMUM  ENTROPY  ACCELERATION  ENERGY  DENSITY 
SPECTRUM  FOR  N4S°E  COMPONENT  OF  STATION  8;  BEAR  VALLEY 
EARTHQUAKE  OF  JUNE  22,  1973  (HIGH  PASS  FILTER  CORNER 

FREQUENCY  AT  0.7  Hz) 
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FREQUENCY  C HZ.  ) 
STATION  8 COMP  N45E 

FIGURE  1-5 


NORMALIZED  MAXIMUM  ENTROPY  VELOCITY  ENERGY  DENSITY 
SPECTRUM  FOR  N45°E  COMPONENT  OF  STATION  8;  BEAR 
VALLEY  EARTHQUAKE  OF  JUNE  22,  1973  (HIGH  PASS 
FILTER  CORNER  FREQUENCY  AT  0.7  Hz) 


STATION  8 

FIGURE  1-6 
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FIGURE  1-7 

NORMALIZED  MAXIMUM  ENTROPY  DISPLACEMENT  ENERGY  DENSITY 
^ SPECTRUM  FOR  N45°E  COMPONENT  OF  STATION  8;  BEAR  VALLEY 
EARTHQUAKE  OF  JUNE  22,  1973  (UNFILTERED  SPECTRUM) 
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Finally,  in  Table  1-2  the  source  parameters  of  the  Parkfield 
earthquake  and  the  Bear  Valley  event  of  June  22,  1973  are  given.  The  most 
important  feature  is  the  maximum  entropy  estimate  of  as  compared  to 

those  given  in  Table  1-1.  In  all  cases,  the  estimate  is  lower,  probably  re- 
sulting from  the  decreased  leakage  from  zero  frequency. 

D.  MODIFIED  MAXIMUM  ENTROPY  SPECTRAL  ANALYSIS 

With  either  maximum  entropy  or  DFT  spectral  estimation,  the 
question  of  how  best  to  incorporate  instrument  response  corrections  for  obtain- 
ing displacement  spectra  from  acceleration  input  data  remains.  Two  standard 
approaches  to  this  problem  are  twice  integrating  the  acceleration  data  in  the  time 
domain  to  yield  a displacement  waveform,  which  is  then  Fourier  analyzed  or 

to  perform  the  integration  in  the  frequency  domain  by  division  of  the  accelera- 

2 

tion  spectra  by  the  function  (2xrf)  . These  methods  are  the  ,re tically  analogus, 
but  division  in  the  frequency  domain  tends  to  be  superior  due  to  accumulating 
errors  which  exist  with  numerical  integration  in  the  time  domain. 

However,  both  method?  cause  errors  in  the  estimation  of  the 

-2 

low  frequency  level,  for  the  function  (2.7Ti)  tends  to  infinity  as  frequency 
goes  to  zero.  Then,  if  the  acceleration  spectral  amplitudes  do  not  tend  towards 
zero  in  a manner  which  would  yield  a limit  when  the  division  is  performed,  the 
low  frequency  level  of  the  displacement  spectrum  will  also  tend  towards  infinity. 
Since  the  behavior  of  the  low  frequency  level  of  the  displacement  spectra  is  im- 
portant in  estimation  of  the  seismic  moment  using  Brune's  model,  and  in  the 
evaluation  of  various  source  models,  such  as  Brune’s  model  and  Archambeau's 
model,  it  is  important  to  eliminate  this  effect. 

A new  technique  of  spectral  analysis  which  should  eliminate  this 
problem  has  recently  been  proposed  by  Robert  Sax  of  Texas  Instruments  Incor- 
porated (Sax,  personal  communications,  1975).  This  method  is  essentially  a 
modification  of  the  maximum  entropy  spectral  estimator  in  which  an  N-point 
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TABLE  1-2 

MAXIMUM  ENTROPY  SOURCE  PARAMETERS 


PARKFIELD 

f (Hz) 
o 

n 

o 

STATION  2 

1.7 

0.  794 

STATION  5 

2.  7 

(6.3  P-Wave?) 

0.  2587 

STATION  8 

3.  3 

0. 1123 

TEMBLOR 

2,7 

(6.  2 P - Wave  ? ) 

0.  190 

I. 59  X 10 


1.  03  X 10' 


1. 22  X 10‘ 


BEAR  VALLEY 
STATION  2 

STATION  7 

STATION  8 


f (S) 
o 

Vp) 

n 

o 

M 

o 

3.  8 

6.  6 

0. 0083 

21 

9.  5 X 10 

2.  8 

6.  8 

0.  025 

22 

3.  1 X 10 

2.  4 

6.4 

0.  0502 

, 22 
6.  0 X 10 
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prediction  filter  is  designed  on  the  signal  autocorrelation  function.  Though 
the  mathematics  of  the  method  does  not  require  it,  the  solution  filter  is  pre- 
sently restricted  such  that  within  the  N-points  the  autocorrelation  function  is 
modeled  only  by  poles.  To  achieve  this,  the  instrument  response,  (27tf)  , 

must  be  removed  as  it  will  introduce  zeros  in  the  prediction  ftlter.  This  tech- 
nique guarantees  a finite  zero  frequency  amplitude  and  zero  amplitude  at  in- 
finite frequency  for  displacement  spectra,  conditions  based  on  Brune's  earth- 

quake  source  mechanism  model 

As  was  previously  discussed,  the  maximum  entropy  spectral 
estimator  generates  a prediction  filter  based  on  the  data  which  can  then  be 
used  to  evaluate  the  spectrum.  The  prediction  filter  is  essentially  obtained  by 

solution  of  the  matrix  equation. 

Rf  = G 

where  R is  the  autocorrelation  spectra  of  the  signal,  f is  the  prediction 
ter  and  G is  the  autocorrelation  function  of  a spike.  It  is  G which  causes 
the  pre-whitening  of  the  spectra.  In  Sax's  modification  of  the  maximum  entropy 
spectral  estimator,  G is  replaced  by  the  autocorrelation  function  of  the  instru- 
ment response,  with  the  poles  removed,  convolved  with  a spike. 

At  present  the  numerical  methods  used  for  this  technique  are 
under  development.  Successful  results  have  been  obtained  on  theoretical  sig- 
nals. When  applied  to  real  data,  though,  the  calculations  have  been  unstable. 

A more  complete  explanation  of  the  theory  and  application  will  be  given  when 
the  formulation  of  the  method  is  complete,  including  a discussion  of  those  tech- 
niques  which  insure  stability. 


1-22 


E.  CONCLUSIONS 

Spectral  estimates  of  near-field  acceleration  data  were  obtained 
using  both  a discrete  Fourier  transform  (DFT)  and  a maximum  entropy  spectral 
estimator.  For  the  former  estimator,  a scaling  error  had  been  found  in  our 
previous  work.  Correcting  this  error  reduced  our  seismic  moment  estimates 
so  as  to  agree  more  closely  with  those  of  other  investigators.  Since  no  change 
in  spectral  shape  occurred,  the  corner  frequency  estimates  were  still  valid. 

Using  the  maximum  entropy  spectral  estimator,  we  obtained 
general  agreement  on  the  spectral  shape  and  amplitude  level  with  the  DFT  esti- 
mate. As  expected,  the  maximum  entropy  spectra  is  considerably  smoother. 
Because  this  estimator  eliminates  leakage,  a more  accurate  estimate  of  the 
higher  frequency  end  is  obtained,  yielding  better  roll-off  values  and  definition 
of  side  lobes.  This  estimator  was  applied  to  data  from  the  Parkfield  earthquake 
and  the  Bear  Valley  evect  of  June  22,  1973;  in  all  cases,  the  low  frequency  level 
was  lower,  probably  resulting  from  decreased  leakage. 

Finally,  work  was  begun  on  a modified  form  of  the  maximum 
entropy  spectral  estimator,  in  order  to  obtain  a better  estimate  of  the  low  fre- 
quency level.  This  method  essentially  involves  an  N-point  prediction  filter 
designed  on  the  signal  autocorrelation  function.  The  technique  guarantees  a 
finite  zero  frequency  amplitude  and  zero  amplitude  at  infinite  frequency  for  dis- 
placement spectra.  Future  development  will  include  the  numerical  methods  used 
for  this  technique  and  application  to  real  data. 


1-23 


I 

t • 

I II 

SECTION  II 

i 

FAR-FIELD  SOURCE  STUDIES 


A.  INTRODUCTION 

For  the  past  several  months,  the  examination  of  far-fie’d 
spectra  for  source  characteristics  has  been  conducted  along  several  lines  of 
investigation.  Surface  wave  data  from  the  Bear  Valley  earthquake  of  June  22, 
1973  were  used  to  determine  a source  mechanism  solution,  and  this  solution  is 
compared  to  that  obtained  from  near-field  data  in  Subsection  B.  The  far-field 
solution  to  the  central  California  earthquake  of  November  28,  1974  is  also  dis- 
cussed in  this  subsection.  In  Subsection  C,  the  results  of  an  experiment  to 
measure  the  effect  of  multipathing  on  surface  wave  magnitude  measurements 
(M^)  and  the  scatter  of  M^-m^  plots  are  discussed.  Theoretical  higher  mode 
spectra  are  discussed  in  Subsection  D,  with  particular  emphasis  on  depth  effects. 
Finally,  in  Subsection  E,  we  summarize  our  results  and  discuss  future  plans. 


B.  FAR-FIELD  SOURCE  MECHANISM  FOR  TWO  CALIFORNIA  EARTH- 
QUAKES 

Using  spectral  fitting  procedures  previously  described  by  Tsai 
(1972),  Turnbull  et  al.  , (1973),  and  Turnbull  et  al.  , (1974a),  we  attempt  to  obtain 
source  mechanism  solutions  to  the  June  22,  1973  and  November  28,  1974 
California  earthquakes  using  the  available  VLPE  and  array  data.  The  event 
description  of  the  former  is  given  in  Table  II- 1,  with  the  travel  paths  to  the 
available  stations  shown  in  Figure  II- 1.  We  see  from  this  figure  that  all  of  the 
travel  paths  are  continental  and  lie  to  the  east  of  the  event.  The  vertical,  trans- 
verse, and  radial  components  recorded  at  these  stations  are  given  in  Appendix  A. 
Demultipathing  procedures  were  applied  to  the  data  on  the  PDP-15  interactive 
graphics  system  (see  Appendix  B). 
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TABLE  II- 1 

EVENT  DESCRIPTION:  THE  BEAR  VALLEY  EARTHQUAKE 

OF  JUNE  22 , 1973 


Event  I.  D.  : BEV/622/73 

Location 

Magnitude 

Date 

Origin  Time 


3 6° 35.  4 N,  121°1 1. 6W 
mb  = 3.  5 
06/22/73 
01:29:12.  3 


Recording 

Location 

Ay.  imuth 

(km) 

Station 

Latitude 

Longitude 

Source 

OGD 

41.  07N 

74. 62  W 

68.  3 

4026.  6 

ALQ 

34.  94N 

106.  46  W 

93.  5 

1343.  6 

LASA 

46.  69N 

106. 22 W 

43.  2 

1672.  5 

TRAVEL  PATHS  TO  THE  AVAILABLE  STATIONS  FOR  THE 
JUNE  22,  1973  EARTHQUAKE 


In  order  to  generate  theoretical  spectra  to  "fit"  to  the  observed 
spectra,  a layered  half  space  is  required  which  provides  a fair  approximation 
of  the  geologic  structure  in  the  source  region.  Although  surveys  are  presently 
being  undertaken  to  determine  this  structure,  no  definite  estimate  exists  for 
this  region.  Helmberger  (1975),  using  shear  body  waves,  has  determined  a 
rough  approximation  for  the  structure  to  the  east  of  the  event  in  the  source 
region.  This  is  shown  in  Table  11-2  below  and  in  Figure  11-2. 


TABLE  11-2 

BEAR  VALLEY  EARTH  MODEL-NORTHEAST  OF  FAULT 


Layer  Number 

Thickness(km) 

V (km/sec) 
P 

V (km/sec) 

p (gm/cc) 

1 

1.  5 

2.  55 

1.5 

2.  0 

2 

13.  5 

5.71 

3.  3 

2.  6 

3 

2.  0 

5.79 

3.4 

2.7 

4 

2.  0 

6.22 

3.6 

2.  8 

i 

21.0 

6.  58 

3.8 

3.0 

* Normal  Gutenburg-Bullen  for  depths  >40  km. 


The  Love  and  Rayleigh  wave  group  and  phase  velocity  dispersion  of  this  model 
is  shown  in  Figure  11-3.  Its  main  characteristic  is  higher  velocities  at  the 
shorter  periods  than  a normal  Gutenberg-Bullen  model  (Turnbull,  et  al.  , 1974a). 

The  spectral  fit  obtained  using  this  structure  is  shown  in  Figure 
11-4.  We  see  that  the  fit  is  of  "average " quality  at  best;  this  is  due  to  the  rela- 
tively small  azimuthal  spread  of  stations  and  the  gross  source  structure  approx- 
imation. The  source  parameter  distributions  are  shown  in  Figure  11-5,  with 
source  parameter  estimates  given  in  Table  II-3.  Considering  the  quality  of  the 
fit,  the  agreement  between  this  solution  and  that  obtained  using  acceleration  data 
(Turnbull  and  Battis,  1974)  is  remarkable.  The  depth,  dip,  slip,  and  strike  are 
in  close  agreement.  The  largest  discrepancy  occurred  between  the  moments, 
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x Observed 


o Theoretical 


FIGURE  II - 4 


SPECTRAL  FIT  - THE  BEAR  VALLEY  EARTHQUAKE 
OF  JUNE  22,  1973 


TABLE  11-3 


ESTIMATES  OF  SOURCE  PARAMETERS  FOR 
BEAR  VALLEY  EARTHQUAKE  OF  JUNE  22,  1973 


Solutions  for  Event:  BEV/622/73 

A.  Solution  by  Minimum-Residual  Criterion 


Optimal  Solution 


Event  I.  D. 

Depth 
h km 

Dip  Angle 

6° 

Slip  Angle 

X° 

Strike 
N 0°  W 

Moment 
, „25  , 

10  dvne-cm 

BEV/ 622/73 

11.  5 

90 

0 

55,  145 

-3 

0.  743  x 10 

Solution  by  Distribution-of-Minimum-Residual  Criterion 


Event  I,  D. 

Probable  Range 

% Confidence 

Source  Parameter 

h(km) 

10.5  - 14.5 

81 

<5 

90° 

80° 

55 

45 

X 

0° 

100 

6 

55,  145 

75 

with  the  near-fieid  moment  being  an  order  of  magnitude  larger.  At  this  time, 
it  is  difficult  to  say  whether  this  difference  is  a result  of  modeling  inaccuracies 
or  a real  source  property.  Future  studies  will  be  concerned  with  a closer 
examination  of  this  problem. 

The  second  California  earthquake,  that  of  November  28,  1974, 
was  analyzed  using  the  same  procedure.  The  event  description  is  given  in 
Table  11-4,  with  the  travel  paths  to  the  available  stations  shown  in  Figure  II-6. 

The  vertical,  transverse,  and  radial  components  recorded  at  thuse  stations 
are  given  in  Appendix  A.  After  demuitipathing  procedures  were  applied,  we 
attempted  to  fit  spectra  recorded  at  stations  wh  se  travel  patho  lie  to  the  east 
of  the  fault.  This  limitation  was  imposed  because  our  source  structure  was 
only  valid  for  that  side  of  the  fault.  Future  studies  will  use  structures  from 
both  sides  to  take  advantage  of  all  of  the  available  data. 

The  spectral  fit  obtained  for  this  event  is  shown  in  Figure  11-7. 

Again,  we  see  that  the  fit  is  of  "average"  quality  at  best  and  for  the  same  reasons. 
The  source  parameter  distributions  are  shown  in  Figure  11-8,  with  the  source 
parameter  estimates  given  in  Table  II-5.  Comparing  our  results  with  those 
obtained  by  Johnson  (1975)  from  body  wave  data,  the  agreement  is  again  quite 
remarkable  considering  the  quality  of  the  fit.  The  depth  is  within  a kilometer 
of  his  solution,  with  the  fault  nearly  vertical  and  mostly  strike-slip  as  is  Johnson’s. 
The  strike  direction  is  approximately  30  degrees  different  using  the  complementary 
solution.  With  a more  accurate  source  structure,  and  utilization  of  all  the  avail- 
able data,  it  is  hoped  that  a better  fit  can  be  obtained.  This  future  result  will 
then  be  compared  to  more  exact  solutions  obtained  from  body  wave  data. 

C.  THE  EFFECT  OF  DEMULTIPATHING  ON  THE  MINIMIZATION  OF 
THE  SURFACE  WAVE  MAGNITUDE  (Mg)  VARIANCE 

In  previous  work  by  Turnbull  et  al.  , (1974b),  the  Mg-m^  dis- 
criminant was  investigated  with  the  objective  of  reducing  the  scatter  of  the 
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TABLE  II-4 

EVENT  DESCRIPTION:  THE  CENTRAL  CALIFORNIA  EARTHQUAKE 

OF  NOVEMBER  28,  1974 


Event  I.  D.  : CCA/332/74 

Location  : 

Magnitude 

Date 

Origin  Time 


36°54.7N,  121°29.  8W 
mb  = 5.2  (MJ 
11/23/74 
23:01:24.  8 


Recording 

Station 


LAS  A 


NORSAR 


TLO 

ALPA 

KIP 

ZLP 


Location 

Latitude 

Longitude 

46.  69  N 

106.22  W 

60.  844N 

10.  887  W 

39.  86  N 

4.  02  W 

65.  233N 

147. 743  W 

21.42  N 

158.02  W 

16.  50  S 

68. 13  W | 

Azimuth 
F rom 


-20.4 


105.  7' 


126.  3l 


A(km) 


1665. 6 


8341.  3 


9380.4 

3589.4 


3914. 0 


8139.  0 


TRAVEL  PATHS  TO  THE  AVAILABLE  STATIONS  FOR  THE 
NOVEMBER  28,  1974  EARTHQUAKE 


o Theoretical 


HQ,Q  10.0  25.0  *+0. 

T C SEC  D 

FIGURE  II -7 

THE  CENTRAL  CALIFORNIA  EARTHQUAKE 
OF  NOVEMBER  28,  1974 


TABLE  II- 5 


ESTIMATES  OF  SOURCE  PARAMETERS  FOR  THE 
CENTRAL  CALIFORNIA  EARTHQUAKE  OF  NOVEMBER  28, 


1974 


Solutions  for  Event:  CCA/  332/74 

A.  Solution  by  Minimum  - Re  s idual  Criterion 


Optimal  Solution 


Event  I,  D. 

Depth 
h km 

Dip  Angle 

6° 

Slip  Angle 

X° 

Strike 

n<£°  w 

Moment 

in25, 

10  dyne-cn 

CCA/332/74 

8.  50 

70 

-30.  0 

| 

30 

0.  522  x 10'1 

Solution  by  Distribution-of-Minimum-Residual  Criterion 


Event  I.  D. 


Probable  Range 


% Confidence 


Source  Parameter 

h(km) 

6.50  - 10.50 

91 

6 


60  - 90 


82 


earthquake  population.  The  method  of  performing  this  reduction  was  to  mini- 
mize the  effect  of  the  source  mechanism  radiation  pattern  and  to  obtain  the 
value  at  the  same  period  for  all  events  and  all  components. 

A set  of  66  Eurasian  earthquakes  recorded  at  three  or  more 
V LPE  stations  during  1972  and  1973  was  selected  for  this  investigation. 

Surface  wave  magnitude  (M  ) measurements  were  first  made  manually  on  these 
records  and  the  results  averaged  for  each  event.  By  averaging,  it  was  honed 
to  reduce  radiation  pattern  effects.  Using  calculated  from  the  vertical 
component  Rayleigh  wave  and  averaged  at  three  or  more  stations,  a standard 
deviation  of  (7  = 0.  354  was  obtained  (ail  values  of  <7  given  here  were  computed 
using  a fitting  technique  which  treats  both  variables  as  independent,  and  mini- 
mizes the  perpendicular  distance  from  the  fitting  line  to  the  data  points).  This 
is  almost  identical  to  that  obtained  by  Lambert,  et  al.  , (1974)  for  from  single 
stations.  Using  calculated  from  both  the  Rayleigh  and  the  Love  waves,  a 
standard  deviation  of  <7=0.  368  was  obtained,  which  is  virtually  the  same. 

Using  the  same  data  set,  Mg  measurements  were  made  by  using 
a spectral  estimate  at  20  seconds  period.  The  vertical  Rayleigh  data  obtained 
in  this  manner  yielded  a standard  deviation  of  <7  = 0.  369.  Using  an  obtained 
from  Rayleigh  and  Love  waves  yielded  a standard  deviation  of  cr  = 0.  363. 

Thus,  the  two  methods  used  previously  (M^  averaging  and  com- 
puting M from  spectral  estimates)  produced  no  improvement  over  the  standard 
deviations  computed  from  single  - station  data. 

Another  possible  way  to  reduce  the  scatter  about  the  mean  value 
lme  could  be  to  demultipath  the  data,  accounting  for  path  effects.  The  method 
used  ;or  demultipathing,  incorporating  a series  of  narrow-band  filters,  is  out- 
lined in  Appendix  B.  Before  processing  of  the  data  was  initiated,  the  66  events 
used  previously  were  reviewed  visually.  Seven  events  were  rejected  due  to 


either  suspicion  of  being  mixed  or  to  the  presence  of  malfunction  (spikes,  etc.  ) 
which  would  interfere  with  the  narrow  bandp.  ss  filtering.  The  net  result  was 
that  the  data  base  was  reduced  to  59  events  consisting  of  221  station-events  (a 
station-event  is  one  event  as  recorded  at  one  station).  The  parameters  of 
these  events  are  listed  in  Table  11-6.  Due  to  time  constraints,  only  the  vertical 
component  of  the  Rayleigh  wave  was  processed  for  each  station  event. 

The  great-circle  travel  paths  from  event  location  to  the  VLPE 
stations  are  shown  in  Figure  11-9  (the  key  to  this  figure  is  given  in  Table  11-7). 

The  event  locations  are  placed  in  the  center  of  regions  which  represent  the  events 
which  were  analyzed.  Only  those  stations  which  recorded  one  or  more  events 
from  a given  region  were  used  in  constructing  the  travel  path  between  that  region 
and  station. 

Since  the  size  of  the  data  base  was  changed,  it  was  necessary  to 

recompute  the  standard  deviations  for  the  previously  reported  three  or  more 

stations  average  M and  the  spectral  estimate  M . We  shall  use  as  a reference 

s s 

the  standard  deviation  for  the  case  of  measured  on  the  vertical  component  of 

the  Rayleigh  wave  and  averaged  over  three  or  more  station  values.  (All 

M data  referred  to  below  was  measured  at  a period  of  20  seconds.  ) The 
s 

following  cases  summarize  the  results  of  this  study. 

• Figure  II- 1 0 illustrates  the  reference  data  set.  The  standard 
deviation  for  this  set  is  a = 0.  283. 

• Figure  11-11  illustrates  the  case  where  was  computed  from 
Love  and  Rayleigh  wave  average  values.  The  standard  de- 
viation of  this  set  is  O'  = 0.285,  which  is  essentially  the  same  as 
the  reference  value  of  O . 

• Figure  11-12  illustrates  the  case  where  M was  computed  from 
spectral  estimates  made  from  the  vertical  Rayleigh  wave.  The 
standard  deviation  of  this  set  is  O = 0.291»  which  is  slightly 
greater  than  the  reference  value  of  O . 
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TABLE  II- 6 

DEMULT1PATHED  EVENT  INFORMATION 
(PAGE  I OF  2) 


o r*  o » f v o 

^ ^ « * * * %%% 

vT  h C T h ^ U'  »C  h h \fl  ^ f**  h h 'O  vC  O'  O1  vC'  ff»  o C.0  OC  'JO  or  <D 

(-*%****%*%%%**%**•.*******  * k * k fc  • 

.c  nuo  j i ^ n vf  3jic>C3avf»0'r3tfvff'air»ir-i«  vc  vo  so 

1<  * » ‘ » * » * • ♦ ‘ ‘ ‘ ‘ ‘ ‘ ‘ ‘ ‘ ‘ ‘ ‘ ‘ * 

10  Ol  “ N (\  (N  (\  N f M f'l  M "1  n O')  f'1  fo  '’I  IO  *■  ^ «”  -•  -1  a »~  CM  fM  (N  DO 

Iflfrft.Cl.O'Oill'BiiJfjfuC'Oir-fcO'OfCjCLtOiOjfjOiOiOjOiPjCidiO.fli 


C;  n tr  ? z w -T«5  I a n ?.  or  o s r.  o ^ in  io 

pq  «o  ;j  i-  (N  | a ro  r-  lTi  a IT  ' 'l  N 

Q m a 


r rg  ;r  :?:  a ;?:  o <N  r-  a 
x>  a n:  cm  *“ 


eoc'trmf'jjO'C'ir'-irNh^fN^m'.'K'JsOMnrcNiridir^ 



ifi  rj  a in  ^ ? u‘  r a it  ^ l'>  it  a m ir  in  g in  in  a in  in  uuT'  in  a in  in  ir, 

to  "1  IT  C t-  f-  -J  W N M O ;l  i£  O O'  >"  'Cl  N f • (M  UU.O  IC  IT*  O ^ IT  O ID  '■*"> 

KOmf'.^CN'rC-aM^NI'I^O'ONf'O'fSeM'IT'Cr" 

o m fN  ■&  a cn  r\i  o (m  m a o c;  m x a a r*  ip  ^ m n in  it  o o o d ro  i t fM 


^.0'Ov00'ar'i\00'0'oaa<v>r*av£ia(N*“»— ooaoof'rr>i  f'*o  x 

, , ! , | ) « | | ( | | | I | | t I • I • ) I * • ■ • I I • < 

(-  o ifiir:  <-  hn  ic  c x c 'i  o c o a a i n o <•  a to  it.  ^ a f f»i  m 'ni'1 

^ n i"  c it  n n s o a rj  n <-  in  a a n ir,  a in  n a if.  it  a i"  a n n N 

► J 

fNOXroxr  a r tr  n o cm  c r ir  n re  o »-  m o «-  co  r~  a cr  x u ^ a 

iri  r>i  ir>  bi  m a a is  a n it  r-  in  a a i/i  r-  o bj  'ni  l' t-  ■■•  o >-  n a r 

E-  no  r-  c a r-  o x x O'  o m x r~  x a a a rg  r-  a <N  at  x r > rg  o rg  o m i— 

i n ip  n :>i  c o oiui  t-  rs  n o n I'm-  i’  in  »■  r a iT  bj  a in  t-  in  a n a 

......: * 

ip,  nor^xrooi'i'- (Nr^aroT—  loo'o  r^iooco'rgxvCvoaP'. 

OT-fNT-r-(NiT-orgooT-o<'iT-i-cccoooo<-oo’-Or-o 

(N  (N  (N  IN  M a (N  (N  Ol  N N f J IN  g H IN  r-l  N H IN  IN  IN  M N N r>l  IM  N N 

t-~  f“~  t"  r-  t"  i"'  r-  t'-  t'-  t"-  t~~  r~~  r~  T'  r--  t'  r—  f' r~~  r*  r" 

W \ \ v V WWWWWWS^NWW^VNNVV 
Pi  o in  . ^ min  m o bi  n m y it  <-  a o1  n id  h co  o >■  h o t-  o vo  r~  x 

«T,»-T-«-'-T-(NirirT'ooc'OC}'—  «*-  rg*-  «“  o o <4  »-  raoco 

r,  \ \ \ \ v \ v \ \ \ v vvwwwwvv^vvvvv 

r-  r-  f r-  f-  f »-»-  rg  f j r\i  rg  rg  c.  rg  rj  ri  ro  m m3  <a  t*-  r-  cc  oc  or  r-  »-  »- 
OOOOOOOOCtOOOCtOOOOOOOOOC:-1  OX) 


tr-  b mo  h f-  m f- 

a a in  in  tr  o <-  in 

OCOOOt-r-T- 

I I I I I I I I 

to  '.)  [O  ID  7 < H fO 

U an  to  it.  :*  h n;  u 

w d*  h o h n ’j  ti 

f/l  s;  m i-  «r  n to  nr. 

VO  f.  U f-4  E-  t-l  X.  IT 

I I I I I I I I 


x x v x x 
►a  »-i  *..t  »-♦  >4 


xxx 
O |J  'J 


i in  ro  f'  in  ic 

n ID  if  (D  tc 

r*  t-  t-  fr- 
ill I I 

E-<  pr  X .4  fj 

ro  !(  U H H 

: to  w ?.  X rr 

i h »t  ><  C C 

< E-*  £-•  !»t  S£ 

I I I I I 

: X X X X X' 

! 4 f rl  rl  O 


no  x m-  t~~  in 

a iTi  n in  o 

a ld  io  M> 
l i i i i 

PI  >4  X . 4 O 

(-  H U SD  C1 

IT  ;V.  r H !D 

o;  rr.  *z  < o 

^ x x r 
i i i i i 

O X X X X 

O iJ  4 O rl 


MD  ID  O Dl  X 
ID  C f O ' l)  0 
r~  t~  X -o  or  x 
i i i i i i 

, 4 O S'r4.a;S~ 
V.  IT  rt  H > rt 
H 7.  ID  K H !1 
•X  O PI  O «C  M 

x r o v.  f--  h 

i i i i i i 

X X X X X X 

► 1 rJ  tj  |4  ►.  rJ 


id  mj  niNni-MtoiinDiD  id  a r u)  a «■■  oniinir  f n o r id 
v ">  oi  n ■ f,  a vo  t*-  a x x cc  a ■ «—  eg  m ro  a an  a-  a it  f-  c o r~  ld  *-*-»*■ 
pooco’-l  OCOOOCOorrrrDlDtrinaiD  vTXvor  r*'  C" 
to  C3  O O O O O O O O O O O O O O O O O O O O O O O O O O C O 


)7iiU  LX-TAI  VA-9U0  '!V21/72 


TABLE  II- 6 

DEMULTFATHED  EVENT  INFORMATION 
(PAGE  2 OF  2) 
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TABLE  II - 7 
KEY  TO  FIGURE  II- 9 


LIST  OF  V LPE  STATIONS 


Sf  ation 
Number 


Name 


Charters  Towers,  Australia 
Cniang  Mai,  7 hailand 
Fairbanks,  Alaska 
Toledo,  Spain 
Eilat,  Israel 
Kongsberg,  Norway 
OgdensDurg,  New  Jersey 
Kipapa,  Hawaii 
Albuquerque.,  New  Mexico 
Matsushiro,  Japan 


Station  Coordinates 

Latitude 

20.  09S 

146.  26E 

18.  79N 

98.  98E 

64.  90N 

148.  01  W 

39.  86N 

4.  02  W 

29.  55N 

34.  96E 

59.  65N 

9.  59E 

41. 07N 

74. 62  W 

21.42N 

158. 02 W 

34.  94 N 

106.  46  W 

36.  54N 

138.  21E 

LIST  OF  SEISMIC  REGIONS 


Station 

Designator 


Name 


Greece  - Turkey 
Caspian  Sea 
Southern  Ira,* 

Central  Asia 
Taiwan 
Kurile  Islands 
Southern  Kamchatka 
Northern  Kamchatka 


Coordinates  of  Region  Center 

Latitude  ! 

Longitude 

37.  IN 

24.  7E 

37.  7N 

51. 2E 

28.  6N 

52.  6E 

35.  6N 

86.  2 E 

22.  ON 

121. 4E 

44.  2N 

148.  8E 

51.  ON 

157.  3E 

55.  3N 

162.  0E 

Figure  11-13  illustrates  the  case  where  Mg  was  computed  from 
spectral  estimates  made  from  both  the  l^ve  and  vertical  Rayleigh 
waves.  The  standard  deviation  of  this  set  i cr  = 0.  276,  which  is 
slightly  lower  than  the  reference  value  of  CT  . 


Figure  11-14  illustrates  the  case  where  Mg  values  were  computed 


by  the  method  of  demultipathing  as  described  in  Appendix  B.  These 
Ms  values  are  averages  of  all  available  station  Mg  values.  The 
standard  deviation  of  this  set  is  a = 0.291,  which  is  slightly  greater 


than  the  reference  value  of  <7. 


Figure  11-15  illustrates  the  case  where  an  effort  was  made  to  re- 
move the  high  and  low  station  Mg  values.  For  each  event,  the 
station  M values  derived  from  the  demultipathing  method  were 
averaged  and  a standard  deviation  computed.  Those  station  Mg 
values  which  were  greater  than  one  standard  deviation  from  the 
average  in  absolute  value  were  rejected.  The  average  of  the  re- 
maining station  Ms  values  was  then  computed.  The  standard  de- 
viation of  the  resulting  data  set  is  a = 0.27  6,  which  is  slightly 
lower  than  the  reference  values  of  cr. 


In  addition  to  the  above,  standard  deviations  for  the  corresponding 
data  sets  at  30  second  and  40  second  periods  were  computed.  These  showed 
essentially  no  change  from  that  of  the  20  second  data. 


The  above  discribed  data  sets  used  mb  values  as  reported  by  PDE. 
Standard  deviations  for  the  data  sets  were  also  computed  using  NORSAR  and 

LASA  m values.  These  standard  deviations  were  markedly  higher  than  those 
b 

using  PDE  m^  data. 


Therefore,  we  have 

the  scatter  in  M^-m^  pl°ts  as  ^eSC 
that,  for  the  case  where  PDE  mfa  v 


seen  that  the  various  methods  of  reducing 
ribed  above  have  not  succeeded.  This  suggests 
alues  are  used  and  all  mixed  events  have  been 
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removed,  a standard  deviation  of  approximately  <7  = 0.28  is  the  minimum  for 

the  data  points  of  an  M -m,  plot  and  cannot  be  reduced  by  attempting  to  reduce 
r s b 

the  variance  in  M . If  any  reduction  in  the  scatter  of  these  plots  is  to  be  achieved, 
s 

it  must  be  accomplished  through  reduction  of  the  variance  of  m^  It  is  suggested 
that  this  be  attempted  in  the  future  by  means  of  maximum  likelihood  estimation 
(Ringdal,  1979).  The  close  agreement  of  the  lines  fitted  tc  the  M^-rr^  sets  pro- 
duced by  the  above  methods  and  the  similarity  of  the  associated  standard  devia- 
tions indicates  that  the  demultipathing  method  is  a good  way  to  determine 
values.  It  has  an  advantage  over  the  presently  used  manual  measurement  method 
in  that  it  is  not  always  possible  to  measure  M^  at  all  of  the  desired  periods  manu- 
ally, whereas  these  measurements  can  always  be  made  by  the  demultipathing 
technique . 
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D EXAMINATION  OF  THEORETICAL  HIGHER  MODE  SPECTRA 

* o 

Using  the  Harkrider  (1964)  formulation  of  a double  couple  source 
in  a layered  half  space,  theoretical  first  higher  mode  Rayleigh  and  Love  wave 
spectra  were  generated.  The  objective  of  this  study  is  to  use  higher  mode  informa- 
tion to  determine  depth  and  structure  in  the  source  region,  especially  for  shallow 
events.  We  can  summarize  the  properties  of  this  spectra  and  its  use  as  follows: 

• Because  the  higher  mode  energy  is  greater  at  shorter  periods, 

its  spe:tra  offers  higher  resolution  than  fundamental  mode  sur- 
face wave  data  for  excitation  differences  between  shallow,  thin 
layers.  As  an  example,  the  variation  of  this  spectra  with  depth 
is  shown  in  Figures  II-16a  and  16b  for  a vertical  strike-slip  fault. 
While  the  spectral  shape  is  similar  to  that  of  the  fundamental  mode 
Rayleigh  wave  of  longer  periods,  the  higher  mode  Love  wave 
spectra  contains  a great  deal  more  character  than  its  fundamental 
counterpart  (Turnbull,  et  al.  , 1974a).  Therefore,  for  the  first 
higher  mode,  both  Rayleigh  and  Love  waves  are  diagnostic  with 
source  depth. 
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FIGURE  II- 1 6a 

FIRST  HIGHER  MODE  RAYLEIGH  WAVE  SPECTRA 
FOR  A VERTICAL  STRIKE-SLIP  FAULT 
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FIGURE  II- 1 6b 


FIRST  HIGHER  MODE  LOVE  WAVE  SPECTRA 
FOR  A VERTICAL  STRIKE-SLIP  FAULT 
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• There  are  several  difficulties  involved  in  the  use  of  higher  mode 

data.  Because  we  are  interested  in  shorter  period  -nergy  (less 
than  15  seconds),  attenuation  and  microseismic  noise  become  a 
problem.  Therefore,  one  is  usually  restricted  to  using  data 
from  very  large  earthquakes  which  generate  long  period  higher 
mode  energy.  In  the  past,  we  have  also  been  restricted  by  the 
instrument  response,  with  very  little  energy  recorded  below  10 
seconds  period.  Some  of  these  difficulties  will  be  eliminated 
though,  when  SRO,  ILPA,  and  certain  experimental  LRSM  data 
becomes  available.  These  installations  utilize  new  broadband 
instrumentation,  with  a fiat  response  to  below  3 seconds  period. 
Employing  filters  for  the  microseismic  noise  and  first  zone 
events  (<20°),  higher  mode  data  should  be  available  over  a range 
of  magnitudes. 

Our  software  for  fitting  fundamental  mode  surface  wave  spectra 
is  presently  being  modified  to  handle  higher  mode  data.  When  this  is  complete, 
and  the  broadband  aata  is  available,  we  will  attempt  to  determine  the  source 
mechanism  and  source  region  structure  using  higher  mode  data. 


E.  CONCLUSION 

Several  tasks  have  been  undertaken  in  our  study  of  far-field 
spectra  for  source  characteristics.  Using  spectral  fitting  procedures  on  funda- 
mental mode  surface  wave  data  from  two  central  California  events,  generally 
close  agreement  was  obtained  with  source  mechanism  solutions  from  body  waves. 
The  spectral  fits  were  of  "average"  quality  at  best,  which  was  due  to  the  rela- 
tively small  azimuthal  spread  of  stations  and  the  gross  approximation  of  source 
region  structure.  For  the  Bear  Valley  earthquake  of  June  22,  1973,  the  seismic 
moment  obtained  from  the  surface  wave  data  was  an  order  of  magnitude  smaller 
than  that  obtained  from  acceleration  data.  It  is  difficult  to  say  whether  this 
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difference  is  a result  of  modeling  inaccuracies  or  a real  source  property.  For 
the  central  California  earthquake  of  November  28,  1974,  oar  solution  was  limited 
because  of  the  availability  cf  earth  structure  information.  Future  studies  to 
these  two  events  will  involve  implementation  of  more  accurate  earth  structures 
for  a definiti\e  comparison  with  bodywave  data. 

In  an  attempt  to  reduce  the  scatter  in  M -m,  plots  for  an  earth- 

s b 

quake  population,  several  attempts  were  made  to  reduce  the  variance  in  M . 

Using  a demultipathing  procedure  as  implemented  for  the  radiation  pattern  with 

four  station  averages,  produced  no  noticeable  effect.  If  any  reduction  in  the 

scatter  of  these  plots  is  to  be  achieved,  the  variance  of  m must  be  reduced,, 

b 

Two  byproducts  cf  this  investigation  were  the  collection  of  enough  source- 
station  travel  path  data  for  a group  velocity  study  of  Eurasia,  and  an  accurate 
method  to  determine  using  the  demultipathing  procedure. 

Finally,  theoretical  higher  mode  surface  wave  spectra  ha\e  been 
generated  for  a double  couple  source  in  a layered  half  space.  Doth  the  Rayleigh 
and  Love  wave  higher  mode  spectra  vary  shape  to  a higher  degree  as  a function 
of  source  depth;  they  are  quite  sensitive  to  source  depth  and  source  region  struc- 
ture. Upon  modification  of  our  spectral  fitting  programs,  and  imminent  avail- 
ability of  broadband  data,  we  should  obtain  better  depth  resolution  for  shallow 
events. 
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APPENDIX  A 

LONG-PERIOD  SEISMOGRAMS  FROM  TWO  CALIFORNIA 

EARTHQUAKES 


In  this  appendix,  the  long-period  seiomograms  used  in  the  source 
mechanism  determination  of  the  22  June,  1973  and  the  28  November,  1974 
central  California  earthquakes  are  displayed.  Each  of  the  figures  represent 
copies  of  the  signal  displayed  on  the  PDP-15  interactive  graphics  system 
(Turnbull,  1974o).  The  lines  at  the  top  of  each  figure  give  the  event  infor- 
mation and  then  the  station  identification.  At  the  bottom  of  each  figure,  the 
amplitude  of  the  trace  is  given  in  milli-microns  (it  is  denoted  by  NM  on  the 
figure).  The  interval  between  scale  marks  on  the  baseline  represents  100 
seconds.  The  lines  denoting  S,  LQ,  LR,  andLREND  indicate  arrival  times 
generated  from  standard  travel  time  curves,  and  are  meant  as  an  aid  to  the 
analyst  rather  than  identification  of  the  particular  phases.  The  seismograms 
displayed  are  as  follows: 

• Bear  Valley  Earthquake  - Junef22,  1973 

Figures  A-l-a,  b,  c:  Vertical,  transverse,  and  radial  LP 
components  as  recorded  at  LASA 
Figures  A-2-a,  b,  c:  The  three  LP  components  as  recorded 
at  OGD 

Figures  A-3-a,  b,  c:  The  three  LP  component!  as  recorded 
at  ALQ 

• Central  California  Earthquake  - November  28,  1974 
Figure!  A-4-a,  b,  c:  Vertical,  tran!verse,  and  radial  LP 

components  as  recorded  at  LASA 
Figures  A-5-a,  b,  c:  The  three  LP  components  as  recorded 


A-l 


at  NORSAR 

Figures  A-6-a,  b,  c:  The  three  LP  components  as  recorded 

at  AFP  A 

Figures  A-7-a,  b,  c:  The  three  LP  components  as  recorded 
at  TLO 

Figures  A-8-a,  b,  c:  The  three  LP  components  as  recorded 

at  KIP 

Figures  A-9-a,  b,c:  The  three  LP  components  as  recorded 

at  Z LP . 
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THE  LONG-PERIOD  TRANSVERSE  COMPONENT  OF  THE 
NOVEMBER  28,  1974  CENTRAL,  CALIFORNIA  EARTHQUAKE 

AS  RECORDED  AT  NORSAR 
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FIGURE  B-l 

DISPLAY  OF  A SINKIANG  EARTHQUAKE  AS  RECORDED 
ON  THE  VERTICAL  COMPONENT  AT  KON 
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FIGURE  B -2 -a 

THE  EFFECT  OF  NARROW  BANDPASS  FILTERS  ON  THE  RAYLEIGH  WAVE 
ON  THE  SINKIANG  EVENT:  10  TO  35  SECONDS  PERIOD 
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FIGURE  B-2-b 

THE  EFFECT  OF  NARROW  BANDPASS  FILTERS  ON  THE  RAYLEIGH  WAVE 
ON  THE  SINKIANG  EVENT:  40  TO  65  SECONDS  PERIOD 
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FIGURE  B-3-a 

THE  ENVELOPE  FUNCTION  OF  THE  BANDPASS  RESULTS  WITH  PEAK 
CHOICES:  10  TO  35  SECONDS  PERIOD 
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FIGURE  B-3-b 


THE  ENVELOPE  FUNCTION  OF  THE  BANDPASS  RESULTS  WITH  PEAK 
CHOICES:  40  TO  65  SECONDS  PERIOD 
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,r)  KEY  ttOROS  fConfinuo  nr.  reverse  side  if  necessary  at  d identity  by  bfork  nwrbei  i 


Seismology 
Near-Field  Spectra 
Maximum  Entropy  Spectra 
Demultjpathing 


Far-Ficld  Spectra 
Higher  Mode  Spectra 
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20  AOSTlVX^T  \l  ontlnuc  on  reverse  ~*dc  !f  nacotsary  and  Identify  by  bh,rk  number;  ~ ~ 

Examination  of  the  seismic  source  from  both  the  near-  and  far -fields 
has  been  undertaken.  Using  a discrete  Fourier  transform  ( DFT)  and  a maxi- 
mum entropy  spectral  estimator  on  near -field  acceleration  data  produced  * 
corner  frequency  and  low  frequency  level  estimates.  This  latter  estimator 
eliminates  leakage,  thereby  obtaining  more  accurate  values  for  the  high  fre- 
quency end,  with  better  roll-off  valuer,  and  definition  of  side  lobes.  When  the 
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20.  continued 

maximum  entropv  spectral  estimator  was  applied  to  data,  from 
the  Parkfield  earthquake  and  the  Bear  Valley  event  of  June  22, 
1973,  the  lower  frequency  level  was  lower  than  that  from  the 
DFT  in  all  cases.  a 

Using  spectral  fitting  procedures  on  fundamental  mode 
surface  wave  data  from  two  central  California  events,  generally 
close  agreement  was  obtained  with  source  mechanism  solutions 
from  bodywaves.  For  the  Bear  Valley  earthquake  of  June  22, 
1973,  the  seismic  moment  obtained  from  the  surface  wave  data 
was  an  order  of  magnitude  smaller  than  that  obtained  from  ac- 
celeration data,  rf-irf  an  effort  to  reduce  the  scatter  in  Mg-  mb 
plots  for  an  earthquake  population,  several  attempts  were  made 
to  reduce  the  variance  in  Ms.  Both  demultipathing  and  correct- 
ing for  tiie  radiation  pattern  had  no  effect.  It  was  concluded  that 
the  variance  of  m^  is  the  controlling  factor.  Finally,  theoreti- 
cal higher  mode  surface  wave  spectra  was  generated  for  a dou- 
ble couple  source  in  a layered  half  space.  Both  the  Rayleigh 
and  Love  wave  higher  mode  spectra  were  found  to  vary  shape 
dramatically  as  a function  of  source  depth. 
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